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ABSTRACT 

We use the dual description proposed by Seiberg to calculate the pressure in the low 
temperature confined phase of N = 1 supersymmetric QCD using perturbation theory to 
O(g^), where g m is the gauge coupling in the dual theory Combining this result with the 
usual high temperature expansion based on asymptotic freedom, we study how the physics 
in the intermediate temperature regime depends on the relative size of the scale parameters 
in the two descriptions. In particular we explore the possibility of having a temperature 
range where both perturbation expansions are valid. 
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1. Introduction 

A standard method to investigate the phase diagram for a theory at finite temperature 
is to use high-temperature and low-temperature expansions of the free energy: At each 
temperature one calculates the free energy by using the variables appropriate at that 
temperature. In QCD, asymptotic freedom tells us that the high T expansion can be done 
in perturbation theory and that the relevant degrees of freedom are quarks and gluons. 
At low temperature, however, the theory becomes strongly coupled, confinement sets in 
and other variables have to be used. The low T degrees of freedom are the hadrons, and 
one calculates the low T expansion using a phenomenological description of a gas of the 
lightest hadrons, i.e. the pions.[] 

Recently, there has been a break-through in the understanding of the infrared structure 
of certain 4D supersymmetric Yang-Mills theories. The strongest results are for N = 2 
SUSY 0, where the low-energy effective action can be calculated exactly, but there are 
also very interesting results for N = 1 P, |J. (For an introduction and a list of references 
see ||). In particular, Seiberg has argued that there often exists two dual descriptions of 
the same theory, one "electric" and one "magnetic". As in usual electric- magnetic duality 
in the presence of monopoles, strong coupling in one description implies weak coupling 
in the other, and a "fundamental" field in one description can appear as a soliton in the 
other. 

In the models considered by Seiberg, the duality between weak and strong coupling is 
manifested in a rather surprising way: Both descriptions are SUSY gauge theories, with the 
same global SU(Nf) flavour symmetry, but with different gauge symmetry. In the case of 
an SU (N c ) color symmetry in the electric theory, the magnetic gauge group is SU(Nf—N c ). 
Since the beta function is ~ (3N C — Nf), the theories are asymptotically free (UVF) or 
infrared free (IRF) depending on the value of Nf. For Nf < \N C) the electric theory is 
UVF and the magnetic IRF. For Nf > 3N C , the reverse holds, and for \N C < Nf < 3N C , 
they are both UVF. Seiberg's duality conjecture, which is supported by several very strong 
consistency checks, asserts that the two theories are identical in the infrared. Furthermore, 
stability of the ground state imposes the condition Nf > N c , and it turns out the the cases 
Nf = N c and Nf = N c + 1 are special 0. In the following we shall restrict Nf to be larger 
than N c + 1. If the gauge group of the electric theory is SO(N), the magnetic theory has 
gauge group SO(Nf — N c + 4) and we shall consider the intervals N c < Nf < 3(N C — 2)/2 
and N c > 4, or Nf = N c > 5, where again the magnetic theory is IRF. It is of interest 
to study both the SU(N) and the SO(N) case since we expect a phase transition from a 
low temperature confined phase to a high temperature deconfined phase for SO(N). In 
the SU (N) case, no deconfining phase transition is expected since there is matter in the 
fundamental representation that will give a perimeter behaviour to the Wilson loop even 
at zero temperature (or, equivalently, a non zero value to the Polyakov loop). 

The main idea of this paper is to exploit the duality in the range N c + 1 < Nf < |iV c 

1 A good introduction to finite temperature field theory containing a list of references to the original 
work is the book by Kapusta m . 
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(for SU(N)) to calculate the leading contributions to the pressure in the electric theory 
both at high and low T. Since this theory is UVF we have the standard perturbative high 
T expansion, but at low energy we can now invoke duality and use perturbation theory in 
the IRF magnetic theory. 

In the next section we outline the calculation of the pressure both at low and high 
temperature while referring some details to an Appendix. In section 3 we discuss the range 
of validity of the expansions and to what extent perturbation theory can be used to obtain 
information about the finite T phase structure. In particular we discuss the importance 
of the scale parameters that control the electric and magnetic coupling constants in the 
asymptotic regions of large and high temperature respectively. 

2. Calculations 

In superfield notation || [7j the action for SUSY QCD is given by, 

S = ^Jd 4 xd 2 6Tr{W a W a ) + h.c. + J d A xd%{4>4i) + (fafa)] (1) 

where W a = W aA T A is the SU(N C ) gauge field strength and the covariantly chiral matter 
superfields 0? transform according to the N c representation of the gauge group (a = 
1,...,N C is the SU(N C ) index, i = l,..,Nf is the flavour index). The covariantly chiral 
superfields <fr a i transform according to the N c representation. Both and <fi are needed to 
get a non-chiral theory. In terms of ordinary fields the theory contains N g = N% — 1 gluons 

(A A ), and gluinos, (X aA , \ &A ), A — l, ..,N g , and also N c N f quarks (ip aa , 4> a \ and 
squarks (<p a , (p*, ip a , y3* a ), a = 1,...,N C . The flavour indices on quarks and squarks are 
suppressed. For SO(N), the number of gluons is N g = N C (N C — l)/2, and the matter fields 
are in a real representation. 

Greens functions in SUSY theories are usually calculated using supergraph techniques 
H [7|. At finite temperature, bosons and fermions do not enter the calculation on equal 
footing, and it is better to use a component formalism.^] In the Appendix we define our 
conventions, and give the component form of the Lagrangian that give the Feynman rules 
shown in figs. 6-9. 

We can now use the standard machinery of finite temperature field theory to calculate 
the free energy density JF = — P at temperature T. The calculations in a theory with zero 
masses and chemical potentials are very straightforward, as demonstrated by an example 
in the Appendix. 

To O(g ) the pressure is due to the "blackbody" radiation from the massless degrees of 
freedom. For SUSY QCD with gauge group SU(N C ), the bosonic degrees of freedom are 
the 2N g gluons and the 4N c Nf squarks. There is an equal number of fermionic partners 

2 Technically this difference enters through the different boundary conditions for bosonic and fermionic 
fields in the imaginary, or "temperature" direction. This difference "breaks the supersymmetry" in the 
sense that thermal expectaion values do not satisfy the zero temperature supersymmetry Ward identities, 
as follows already from the supersymmetry algebra § . The theory is of course still supersymmetric since 
thermodynamics is determined by the (T = 0) spectrum which is supersymmetric. 
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to these (gluinos and quarks), and thus the "blackbody" part of the pressure is given by 



C = (2N g + w c N f )(i + 8J go 

= (N? + 2N c N f -l)^ . (2) 

To get the corresponding result for SO(N), one substitutes the value for N g given above 
and also replaces Nf by Nf/2 since quarks and squarks are in real representations and thus 
contribute only half as many degrees of fredom as in the SU (N) case. In the dual magnetic 
description the gauge group is replaced by SU (Nf — N c ) but the number of flavours remains 
Nf. Following Seibergs duality hypothesis, we also need Nf x Nf free chiral multiplets, 
which for S77(iV)adds 2Nj bosonic degrees of freedom (and an equal number fermionic 
ones), and thus for the dual "magnetic" description we get 



r>mag 



2[(N f - N c ) 2 - 1)] + 4(N f - N c )Nf + 2Nf 



2 ll+ iy_^ 

8/ 90 



= tf + ^-M)^ . (3) 

For SO(N) we have a similar expression remembering that the gauge group is now SO(Nf — 
N c + 4) and that the meson multiplet is formed from a symmetric combination of two 
complex fields, and thus has 2Nf(Nf + l)/2 components. 

Comparing (2) and (3), we see that if Nf < 3N c /2 the number of "blackbody" degrees of 
freedom in the electric description is greater than the corresponding number in the magnetic 
description. This is a reasonable result: The same condition Nf < 3N c /2 describes the 
the range of Nf for which the magnetic description is infrared free and the electric theory 
is confining. (We remind the reader that in the region 3N c /2 < Nf < 3N C Seiberg has 
argued that we should have a non-abelian Coulomb phase and that the electric theory loses 
asymptotic freedom at the upper of these limits.) We would certainly expect the number 
of massless "hadrons" to be smaller than the number of constituents in a confining theory, 
and this is exactly what we find. A similar calculation for the case of SO(N c ), where the 
limiting case is when Nf = 3(N C — 2)/2, gives the same result: The number of massles 
degrees of freedom in the electric theory is larger than the ones in the dual magnetic theory 
as long as the electric theory is confining and the magnetic is infrared free. 

The 0(g 2 ) contribution to the pressure is obtained by calculating diagrams of the type 
shown in fig. 1, using the rules given in the Appendix. The result is 



p 2 = -N g (C A + 3T f ) ^ , (4) 



where C A is the Casimir operator for the adjoint representation, and Tf depends on the 
matter content. For SU(N), C A = N c and T f = N f , while for SO(N), C A = (N c -2)/2 and 
Tf = Nf/2. The pressure in the magnetic theory is obtained by replacing 7V C by Nf — N c 
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Figure 1: Some of the ^-contributions 

or Nf — N c + 4 for SU (N) and SO(N) respectively. (Note that the region we consider, the 
extra chiral fields are non-interacting, so they contribute only to the blackbody pressure.) 

In theories with massless bosons, the next contribution is 0(g 3 ), rather than 0(g 4 ), 
due to the infrared singularities related to the ring diagrams exemplified in fig. 2. The 
techniques for calculating these diagrams are well known, and essentially amonts to calcu- 
lating the (thermal) electric mass, m eh for the gluons and the thermal mass, m sq , for the 
squarks. These masses are readily obtained from the standard self-energy diagrams, and 
we find, 



1 



m 



(C A + T f )g 2 T 2 



el 



(5) 



m 



sq 



N 9 
12N, 



9 Z T 



This implies the following result for the 0(g 3 ) contribution to the pressure, 



->el 



\ 2 J 9 \12N, 



N \ 3/2 
J " 9 ^ N c T f 



12n 



(6) 



(4^ + ^-) 3 / 2 N g , and again the result in the corresponding magnetic theory is obtained by 
replacing N c by N f - N c , or by N f - N c + 4. 





Figure 2: Ring-diagrams 



5 



To get the total pressure to 0(g 3 ) in the electric theory we add Pq, Pf and Pf and 
replace g 2 by a runing coupling constant at temperature T, 

2 = 87T 2 

9el b e MT/A el ) 1 ' 

and for the magnetic theory we add the corresponding results and replace g 2 by 

4tt 2 



mag 



u mag 

ln(A mas /T) 



For SU(N), b el = 3N C - N f and b mag = N f - 3N c /2, while for SO(N), b el = 3N c /2 - 
Nf/2 — 3 and b mag = 3N c /2 — 3 — Nf. The expression for g 2 el is valid for T ^> A d and 
the one for <7^ ag for T <C A mag . If we define the ratio x = A mag /A e i, there is a possibility 
of using perturbation theory in both the electric and magnetic theory simultaneously if 
1 < T I A e \ < x. This will be discussed further below. Finally we should stress a very 
important point: Because of supersymmetry, the vacuum energy is identically zero, and 
thus also the pressure at zero temperature. This is not the case for non-supersymmetric 
theories like QCD, where there is a non-zero vacuum pressure (related to the quark an 
gluon condensates, and phenomenologically manifested by the so called bag constant) that 
is important for understanding the finite T phase transitions |l| . 

3. Discussion 

We shall first discuss the validity of the high and low T expressions for the pressure 
derived in the previous section. It is known that in QCD, the perturbative expansion of 
the (T = subtracted) free energy is well defined only to 0(g 5 ), since at 0(g 6 ) it becomes 
sensitive to the non-perturbatively generated 0(g 2 ) magnetic mass. For a good discussion 
of the results to 0(g 3 ) see [Q. Recently the 0(g A ) term in pure Yang- Mills theory has been 
calculated by Arnold and Zahi ||. The conclusion, so far, is that perturbation theory is 
reliable only at very high temperatures, and this is likely to be true also for the theories 
we are considering here. Unfortunately, because of finite size corrections, not even the 
most recent lattice simulations of pure YM theory are good enough to really settle how 
well perturbation theory works. For theories with dynamical fermions, like the one studied 
here, we may have to wait long before reliable lattice results are available. Another point 
which is relevant for this discussion is the contribution to the pressure from topologically 
nontrivial gauge configurations. In particular, for QCD the contributions from instantons 
have been calculated by Pisarski and Yaffe (for a review see |]lOl) who showed that the 



instantons become important only at temperatures too low for perturbation theory to be 
valid. It would be interesting to repeat this calculation in the super symmetric case. 

The presence of two scale parameters, A d and A mag , requires some comment. A e[ is 
the usual renormalization group invariant scale of the asymptotically free (electric) gauge 
theory, and just as in ordinary QCD it sets the scale for the onset of nonperturbative 
phenomena. While less familiar, the situation is very similar in an IR free theory. The 
reason that this is not always stressed in the prototype IR free theory, QED, is that the 
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electron mass provides a natural renormalization point. Since a is very small at this scale, 
nonperturbative effects become important only at distances so short that QED must be 
extended to the full electroweak theory. If, however, the electron had been massless, using a 
running coupling constant and a renormalization group invariant scale would be as natural 
in QED as in QCD. Again the renormalization group invariant scale parameter Aqed would 
tell us at what scale the non-perturbative (short distance) physics becomes important. 

In the theories we consider, naturalness implies that there is only one scale connected 
to the onset of nonperturbative phenomena. However, this does not mean that the A 
parameter occuring in the UV (or electric) or IR (or magnetic) description of the theory 
are necessarily the same. The A-parameters depends on the renormalization scheme and 
could well be different. In our case g e \ becomes large as T approaches A e \ from above, and 
9mag becomes large as T approaches A mag from below. Thus there are two possibilities: 1) 
A mag < A ei . In this case the temperature region in between A mag and A d is not even in 
principle accessible in perturbation theory. 2) A mag > A e [. Here there is a possibility of 
simultaneously using perturbation theory in both the electric and the magnetic description, 
but there is of course no guarantee, or even great hope, that the perturbative expansions 
will be good in the overlap region. At present we have no way of determining which of these 
alternatives is correct - this could only be settled by matching a direct nonperturbative 
calcualtion in the electric theory, e.g. using lattice methods, to a corresponding weak 
coupling calculation in the magnetic theory. Any concrete idea for how to perform such a 
calculation would be very desireable. 

Clearly, the second of the above alternatives is the more interesting one, and we shall 
now investigate it a bit further. In fig. 3 we show the SU(N) 0(g°), 0(g 2 ) and 0(g 3 ) 
contributions to the pressure as a function of ln(T/A e j) for x = A mag /A e i = 50. We 
see that even though the ratio x is taken (presumably unrealistically) large perturbation 
theory looks very poor in the overlap region. The reason is the large coefficients in front of 
the logarithms in the 0(g 2 ) and 0(g 3 ) terms. One might have hoped to see a qualitative 
difference between the the SU(N) and SO(N), since a phase transition in expected for the 
latter, but in fact fig. 3 is typical for both gauge groups. 

An obvious question is whether there is a range of parameters N c and Nf where the 
coefficents of the logarithms become small. The answer is yes, by simply noting that for 
Nf very close to N c the number of magnetic gluons is small and thus also the coupling in 
the magnetic theory. This is most clearly seen in the large N c limit where the 0(g 2 ) and 
0(g 3 ) corrections vanish for Nf/N c — > I. In fig. 4 we show the different components of 
the pressure in the SO(N) case, where a phase transition in expected, for N c = Nf = 30 
and x = 20. In this, admittedly rather extreme case, it is not unreasonable to assume that 
both perturbation expansions can be trusted in the region in the proximity of the arrow. 
The shape of the curves are however not compatible with having a phase transition, where 
the pressure must be continuous, and concavity of the free energy density requires that 
discontinuities of the derivative of the pressure, and hence of P/T 4 , must be such that it 
increases with T. From this we can conclude that the value x = 20 cannot be correct, 
so the perturbative analysis have in fact allowed us to put a bound on the intrinsically 
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Figure 3: Pressure versus logaritm of the temperature for iV c = 9, Nf = 12, A e ; = 1 and 
A mag — xA e i = 50 A e ;. The thin solid lines are the blackbody contributions, the dotted 
lines the 0(g 2 ), the dashed lines the 0(g 3 ) and the thick lines the total contributions. 

non-perturbative parameter A mag /A e i. Also note, that a bound on x also translates into 
a bound on T c /A e i, where T c is the transition temperature. Clearly these bounds are not 
very strict, and we have not made any attempts to systematically explore different values 
for x, Nf and A^ c . 

In conclusion, we have explored the possibility of using perturbation theory to study 
both the high and low temperature regime in N=l SUSY theories possessing Seiberg dual- 
ity. First we noted that the number of black-body degrees of freedom in the electric theory 
is larger than in the magnetic, exactly in the region where the former is UV freee and the 
second IR free, which is what is naively expected in a confining theory. We calculated the 
pressure to 0(g 3 ) in both the UV and IR regime, and discussed various possibilities for the 
intermediate temperature regime. Our results are somewhat disappointing, in the sense 
that the dream of finding a regime where both perturbation expansions were applicable 
was not fullfilled. We did, however, manage to glean some nonperturbative information 
from our calculations by studying the large A^ c limit. 

There are several directions in which this work might be extended that could be of 
interest. We have already mentioned that instanton effects could change our resulsts, and 
so could of course higher order terms in perturbation theory. Another interesting question 
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Figure 4: Pressure versus logaritm of the temperature for N c = Nf = 30, A et = 1 and 
A-mag — x ^-ei — 20A e ^. The lines have the same meaning as in fig. 3. 

is to study the deformed theories, either by including mass terms, or by breaking the global 
symmetries by choosing moduli parameters different from zero. In particular one could ask 
whether finite temperature effects would induce a flow in moduli space, and thus making 
some vacuua unstable against bubble formation. 

Acknowledgements: We thank A. Fayyazuddin, and P. Di Vecchia for discussions and 
N. Seiberg for a useful correspondence. 
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Appendix 

A. Feynman rules 

In components, the lagrangian density (1) takes the form, 

C = -\f a ^F a + i\ aA (D\) a A + ir a (D^) aa - ir a (D^ a 

+{D^) + \D^) a + (D^r(D^): 

+iV2g [(<p*T A ^ a )\ A + $«T A <p)\ A + X A ('ijj a T A (p*) + (<fT A ^)\« A ~ 

-^ 2 E[(^M-(^V)] 2 (a.i) 

where the matrices T A , A = l...N g span the fundamental representation of SU(N C ), and 
where we have suppressed contractions of group indices. The various covariant derivatives 
are defined as follows: 



(DX) a A = 


{^)^[d,5 AB + gf ABC A^}\ b 


(P1>)a = 


(aT a [dX + wA A (T A ) b a ]^ ab 


(D$) aa = 






[dX + wA A (T A ) b a ] Vh 




d^ a - ig<p* b (T A y b A A 


(D^r = 


d^ a - *g^ b (T A y b A A 




idX + wA A (T A ) b M 



(A.2) 

Our Minkowski metric is rj^ = diag(l, — 1, — 1, — 1), and the 2 x 2-matrices cr M and <x M 
satisfy the relations 

oi**^ + o***^ = 2 V ^5$ (A.3) 

From the action we read off the Feynman rules given in figs. 6-9. For the gauge group 
SO(N) the rules are identical except that since the matter multiplets are real, there are 
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no dotted quark or squark lines. In the calculations, this only amounts to an overall factor 
of 2 for the relevant diagrams. 

B. Sample Calculation 

We illustrate the calculation of the g 2 contributions by evaluating the diagram in fig. 
6 involving one gluon and two gluinos. 




Figure 5: Diagram contributing to the pressure to 0(g 2 ). 
Using the Feynman rules and dropping a temperature independent piece we find, 

'i{ap)W5 oG \ ( , fABC , ( i(*q)™6 AA \ 



l Ggg = -|ff 2 /_ dpdq JdkU %{PV) 2 (~K W 



2ig 2 N c N g f_ dpdq J + dki-^\ (2ir)%p + q + k), (A.4) 



where +(-) denotes the momentum integral for a boson (fermion), cf. (|A.6| ) below. The 



second equality follows directly using cx-matrix algebra, and taking traces. To see that this 
integral factorizes we use the delta function to write 2pq = k 2 — p 2 — q 2 and obtain 

X Ggg = ig 2 N c N g (A 2 _ - 2A + A_) = i\g 2 N c N g A 2 + , (A.5) 

where, 

f dp 1 f d?p rioo+s \ f 2 \ 

A ^Lv 2 ^ T ^ii (2^Loo +£ ^ o ^t(^f7J (A - 6) 



Evaluating these integrals yields, 



A+ = -2A_ = ~. (A.7) 



Notice the factor i in ( |A.7|) . It enters because we use Minkowski space Feynman rules 



to extract the free energy. A proper translation back to Euclidean space shows that the 



11 



result of a L-loop diagram evaluated with the propagators and vertices in figs. 6-9 and 
integrations as in the above example, has to be multiplied with (— to give the proper 
contribution to the free energy. 

The above example is typical; all our g 2 -integrals factorize in this manner and are 
expressible in terms of A±. For the ring diagrams that yield the g 3 contributions, (fig. 
2), we must calculate the self-energy corrections to the bosonic propagators with zero 
momentum on the external legs. This last condition again makes the calculation easy and 
the result can be expressed in terms of A±. Again, since we work in Minkowski space, some 
care is needed when extracting the thermal masses from the Greens functions, in order to 
get the correct phase. 
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Figure 6: Propagators for gluon, ghost, right quark, left quark, gluino, right squark and 
/e/t squark, respectively. 
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Figure 7: Gluon-Fermion Vertices 
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Figure 8: Gluon-Squark and Squark-Squark Vertices 
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Figure 9: S quark- Fermion Vertices 
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